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Arrays of identical oscillators can display a remarkable spatiotemporal pattern in which phase- 
locked oscillators coexist with drifting ones. Discovered two years ago, such "chimera states" are 
believed to be impossible for locally or globally coupled systems; they are peculiar to the intermediate 
case of nonlocal coupling. Here we present an exact solution for this state, for a ring of phase 
oscillators coupled by a cosine kernel. We show that the stable chimera state bifurcates from a 
spatially modulated drift state, and dies in a saddle-node bifurcation with an unstable chimera. 
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PACS numbers: 05.45.Xt, 89.75.Kd 

In Greek mythology, the chimera was a fire-breathing 
monster having a lion's head, a goat's body, and a ser- 
pent's tail. Today the word refers to anything composed 
of incongruous parts, or anything that seems fantastical. 

This paper is about a mathematical chimera in which 
an array of identical oscillators splits into two domains: 
one coherent andphase-locked, the other incoherent and 
desynchronized [ll y, b| ■ Nothing like this has ever been 
seen for identical oscillators. It cannot be ascribed to a 
supercritical instability of the spatially uniform oscilla- 
tion, because it occurs even if the uniform state is sta- 
ble. Furthermore, it has nothing to do with the par- 
tially locked/partially incoherent states seen in popula- 
tions of non-identical oscillators with distributed frequen- 
cies 0, 1^ . There, the splitting of the population stems 
from the inhomogeneity of the oscillators themselves; the 
desynchronized oscillators are the intrinsically fastest or 
slowest ones. Here, all the oscillators are the same. 

In this Letter we explain where the chimera comes from 
and pinpoint the conditions that allow it to exist. It 
was first noticed by Kuramoto and his colleagues 0, 0, 
1^ while simulating arrays of limit-cycle oscillators with 
nonlocal coupling. As they emphasize, nonlocal coupling 
is less explored than local or global coupling 
Hol Ull 0. It arises in diverse applications 
ranging from Josephson junction arrays ^3 ^^'^ chemical 
oscillators P, |3, H. to the neural networks underlying 
snail shell patterns jl4L and ocular dominance stripes 

MM 

We study the simplest system that supports a chimera 
state: a ring of phase oscillators P, |3| governed by 
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G{x — x') sin [<f){x, t) — (j){x' , t) + a] dx' 



(1) 

Here (f>{x,t) is the phase of the oscillator at position x 
at time t. The space variable x runs from — tt to tt with 
periodic boundary conditions. The frequency ui plays 
no role in the dynamics; one can set w = by redefin- 
ing (j) ^ (f> + Lut without otherwise changing the form 
of Eq. (0). The angle < a < -I is a tunable param- 
eter. The kernel G{x — x') provides nonlocal coupling 



between the oscillators. It is assumed to be even, non- 
negative, decreasing with the separation |x — cc'l along 
the ring, and normalized to have unit integral. Ku- 
ramoto and Battogtokh P, |3l assumed an exponential 
kernel G{x) oc exp(— but instead we will take 
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Acosx) 



(2) 



where < yl < 1. Simulations show that both kernels 
give qualitatively similar results, but the cosine kernel 
allows the model to be solved analytically. 

Figure QJa) shows a snapshot of a chimera state for 
Eq. The oscillators near x — zLtt are locked and 

coherent: they all move with the same instantaneous fre- 
quency and are nearly in phase. Meanwhile, the scattered 
oscillators in the middle of Figure 1 (a) are drifting, both 
relative to each other and relative to the locked oscilla- 
tors. They slow down as they pass the locked pack, which 
is why the dots appear more densely clumped there. 

These simulation results can be explained [II by gener- 
alizing Kuramoto's earlier self-consistency argument for 
globally coupled oscillators Let denote the angu- 

lar frequency of a rotating frame in which the dynamics 
simplify as much as possible, and let 9 — (f> — fit denote 
the phase of an oscillator relative to this frame. Intro- 
duce a complex order parameter i?e*® that depends on 
space and time: 



Then Eq.lQJ becomes 
80 



G{x~x')e''^'^'''''Ux' 



^^uj-n- i?sin W-e + a] 
ot 



(3) 
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By restricting attention to stationary solutions, in which 
R and O depend on space but not on time (a condition 
that also determines il), Kuramoto and Battogtokh (J 
derived a self-consistency equation equivalent to 



R{x) exp [iQ[x) 



G(x-a;')exp [i<d{x')] 



R(x') 
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FIG. 1: (a) Phase pattern for a chimera state. Parame- 
ters: A = 0.995, /? = 0.18, iV = 256 oscillators. Equation 
Q was integrated using a Runge-Kutta method with fixed 
time step dt — 0.025 for 200,000 iterations, starting from 



[x) — a exp 



where r is a uniform random variable 



on [- 



We took a = 6 and b = 30. (b) Local phase co- 



herence R{x), computed from J^J. Locked oscillators satisfy 
R(x) > A. (c) Local average phase B(a;). 



where (3 = ^ — a and A = uj — fl. 

Equation Q is to be solved for three unknowns — the 
real- valued functions R{x) and 0(x) and the real number 
A — in terms of the assumed choices of and the kernel 
G{x). Kuramoto and Battogtokh P, solved ||SJ) nu- 
merically via an iterative scheme in function space, and 
confirmed that the resulting graphs for R{x) and Q{x) 
match those obtained from simulations of Eq.f^l- Fig- 
ures ^b) and^^c) show the graphs of R{x) and 9(x) for 
the parameters used in Figure ^a). 

The resemblance of these curves to cosine waves sug- 
gested to us that Eq.lO might have a closed-form solu- 
tion. It does. Since the right hand side of Q is a convo- 
lution integral, the equation is solvable for any kernel in 
the form of a finite Fourier series; that is what motivated 
the choice of Q. For this case, R{x) and Q{x) can be ob- 
tained explicitly. The resulting expressions, however, still 
contain two unknown coefficients, one real and the other 
complex, that need to be determined self-consistently. 

The solution proceeds as follows. Let 



Expanding G with a trigonometric identity, Eq. jSj) gives 



113, 



cos X ) cos X 



sin a;') since 



a cos X 



(7) 



where the unknown coefficients a and c are given by 

c = e"^{he'^^) (8) 
a = Ae'^ihe'^ cos x') . (9) 



The coefficient of sinx vanishes in ((TJ, if we assume 
R{x') = R{—x') and 6(a;') — 6(— x'), as suggested by 
the simulations. The assumed evenness is self-consistent: 
it implies formulas for R{x) and 0(a;) that indeed possess 
this symmetry. For instance, R{x) satisfies 

i?2 = (i?e*®)(i?e-*®) 

= (c + a cos x) (c* -|- a* cos x) 

= |c^| -|-2Re(ca*)cosx-l- lapcos^x , (10) 

which also explains why the graph in Fig. ^b) resembles 
a cosine. Likewise, 0(x) satisfies 



tan 8(x) 



R{x)smQ{x) Ci +ai COS. X 

i?(x)c0s8(x) Cr + ftrCOSX 



(11) 



where the subscripts denote real and imaginary parts. 

Another simplification is that c can be taken to be 
purely real and non-negative, because of the rotational 
symmetry of the governing equations. In particular, the 
self-consistency equation ^ is left unchanged by any 
rigid rotation 0(x) 8(x) -I- Oq. Thus we are free to 
specify any value of 0(x) at whatever point x we like. 
We choose 0(-|) — 0. Then i?e'® = c + acosx implies 
— c. Since R is real and non-negative, so is c. 
Hence, we take Ci — from now on. 

To close the equations for a and c, we rewrite the av- 
erages in © and in terms of those variables. Using 
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A - ^A2~i?2(x) 



R c + a* cos X 

and inserting H1U|I into (jHJ and we find 



A — (A^ — — 2car cosx — jap cos^ x) 
c + a* cos X 



(12) 



Mx') = ^^^i^HW) (6) g = / A-(A^-c^-2e..^cos.-H^co.^.)^ ^ 



R{x') ' ' U-^t y c-Ka* cosx 

' (13) 

and let angular brackets denote a spatial average: This pair of complex equations is equivalent to four real 

equations for the four real unknowns c, Or, at, and A. 

, 1 n i\ 7 I The solutions, if they exist, are to be expressed as func- 

(/) = / fix )dx . r ^1 / ^ J ^ 

2tt J__ tions of the parameters p and A. 
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FIG. 2: Region in parameter space where a chimera state ex- 
ists. Solid line, boundary determined by numerical solution of 
Eqs. I12L 1131 : dashed line, approximate boundary obtained 
from perturbation theory (see text). A saddle-node bifurca- 
tion of a stable and unstable chimera occurs at the boundary. 



Figure |2] plots the region in parameter space where 
chimera states exist, computed by solving Eqs. ifT^ . (fO|l . 
with a root-finder and numerical continuation. To obtain 
starting guesses for the four unknowns, we integrated (fT|l 
numerically, then fit the resulting R[x) and Q{x) to the 
exact solutions ifTUI) , ifTTl) , to extract the corresponding c, 
ttr, and Ui, and estimated A directly from the collective 
frequency of the locked oscillators. By sweeping f3 at 
fixed A, we found that the chimera state disappeared 
suddenly when it reached the boundary of the region. 

For deeper insight into the chimera state and its bi- 
furcations, we now solve Eqs. ifT^ . ((T^ perturbatively. 
Figure 121 suggests that we should allow (3 and A to tend 
to zero simultaneously: let A = e, p = /3ie, and seek 
solutions of p2(l . H13|l as e — > 0. Numerical continuation 
reveals that the solutions behave as follows: 



A 

c 

ai 



1 + Aie -I- A2e^ 
1 + cie + C2e^ 



(14) 



where terms of O (e'^) have been neglected. Substituting 
this ansatz into H12|l . H13() . we find that Ai = ci is re- 
quired to match terms of O (\/e). Then at leading order, 
the real and imaginary parts of Eqs. 112|l . (|13|l become 
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-Re 
Im 



= -Rc 



\/2 ^cos xV 6 ^ u cos X 



-Im \/2 ^cos x\/ 5 ~ u cos x 



(15) 
(16) 
(17) 
(18) 



where we've defined 5 — A2 — C2 and u — a2r- 

The solutions of these equations can be parametrized 
by (5, as follows. Writing /(w, 5) for the right hand side of 
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FIG. 3: Fraction of oscillators in the chimera state that drift. 
Solid line, stable chimera; dotted line, unstable chimera. 



Eq. H17|l . we compute all the real roots of u = f{u, (5), and 
regard them as functions of 5. Then we sweep through 
all the i5's for which Eq. (|17|l has a solution, and substi- 
tute the associated u{8) into the remaining equations to 
generate values of ci((5), /3i((5), and a2i{5). 

This approach gives a great deal of information 
about the chimera state. For instance. Figure O plots 
the fraction of oscillators that are drifting, /drift = 
^ cos~^((5/m((5)), as a function of the control parame- 
ter /3i((5). There are two branches of solutions. The 
upper branch (which dynamical simulations of Eq. 
show to be stable) bifurcates from a state of pure drift 
at /3i = 0. As /3i increases, drifting oscillators are pro- 
gressively converted into locked ones, eventually reach- 
ing a minimum of about 44% drift at the largest j3i for 
which stable chimeras exist, (/3i)„iax ~ 0.2205. There 
the upper branch collides with the lower (unstable) one, 
which itself emerges from a homoclinic locked state at 
(3i = 0, where the in-phase oscillation of Eq. ^ is lin- 
early neutrally stable. The maximum value of /3i predicts 
that the slope of the stability boundary in Fig. [3 equals 
1/ (/3i)inax ~ 4.535, shown there as a dashed fine. 

Unfortunately, when the variables are plotted versus 
j3i, much of the bifurcation structure is hidden. In par- 
ticular, two crucial events in the genesis of the chimera 
state occur when /drift = 1; /^i = 0, and therefore col- 
lapse onto a single point in Figure|21 These events create 
the x-dependence in the chimera state, first in its local 
coherence R{x) and then in its local average phase Q{x). 
To see how such spatial structure arises, it is best to 
treat 5, not /3i, as the relevant parameter, even though 
it is not a true control parameter (its turning points do 
not signify bifurcations, for example). 

Figure ^ plots u{5) vs. 5 for the roots of Eq. l(T7|l . 
Branches have been coded with different dashing styles to 
indicate that they represent qualitatively different states. 
The zero branch along the 5-axis, shown as a solid line, 
represents a family of spatially uniform drift states where 
both R and are independent of x. Such states occur 
only when /3 = 0. They correspond to the exact (non- 
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FIG. 4: Solutions of Eq. (1171 in {S, u) plane, where 5 — /S.2 — C2 
and u = a2r ■ Inset panels show typical mean drift frequencies 
A(a;) in different regions of solution space. Panels with arrows 
indicate the shape of A(a;) for transitional values of {5,u). 



perturbative) solutions of Eq. ||SJ) given by 8(a;) 
and R{x) = R ^ V2A - 1, with \ 







R = V2A - 1, with i < A < 1. Lin- 
earization shows that this uniform drift state undergoes 
a zero-eigenvalue bifurcation at A = i(2 + e)^^, valid for 
all < e < 1, implying a critical value of (5 = ^ as e — > 0. 

This is the first crucial event. At (5 = |, a spatially 
modulated drift state is born. Now R depends on x. In 
perturbative variables, a root with u = a2r 7^ bifur- 
cates off the zero branch (shown dashed in Figure 0J . 
Meanwhile, a2i = 0, so we still have Q{x) = {) for all x 
(from Eq. (|1I|I and Ci =0). Thus for states on the dashed 
branch, all the oscillators are drifting, and maintain the 
same average phase, but with different amounts of co- 
herence at different values of x. Like the uniform drift 
states, these modulated drift states occur only if /3 = 0. 

The second crucial event occurs when the dashed 
branch intersects the line u = 5. Then several things 
happen. The first locked oscillators are born; a2i and /3i 
become nonzero; Q depends on x; and a stable chimera is 
created. Evaluating the integral in (|17|l for u = 5 shows 
that 5 = 16/97r'^ w 0.18 at the birth of the chimera. 

Another way to distinguish among the various states 
is shown in the insets of Figure 0] For selected values 
of 5 we have plotted the time-averaged frequency A(a;) 
of the oscillator at x, measured relative to the rotating 
frame. For locked oscillators, ^{x) = 0; for drifting 
ones, A{x) — ^A^ — R'^[x) ~ t^J 5 — ucosx, to lead- 
ing order in e. Starting from the origin of Figure ^ 
and moving counterclockwise around the kidney bean, 
the corresponding graph of /S.{x) is zero for the homo- 
clinic locked state; flat for uniform drift states; modu- 
lated and positive for modulated drift states; and par- 
tially zero/partially nonzero for chimera states, with the 
fraction of drifting oscillators decreasing steadily as we 



circulate back toward the origin. 

Although we have focused on the chimera state in 
Eq. 1^, it also arises in other spatially extended systems. 
Indeed, it was first seen in simulations of the complex 
Ginzburg-Landau equation with nonlocal coupling Q, 0] ■ 
That equation in turn can be derived from a wide class 
of reaction-diffusion equations, under particular assump- 
tions on the local kinetics and diffusion strength that 
render the effective coupling nonlocal [iL y, |l2| . 

In two dimensions, the coexistence of locked and 
drifting oscillators manifests itself as an unprecedented, 
bizarre kind of spiral wave — one without a phase singu- 
larity at its center 0, |^ . Perhaps our analysis of its one- 
dimensional counterpart can be extended to shed light 
on this remarkable new mechanism of pattern formation. 
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